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SIMPLIFIED CALCULATION OF TRANSITION MATRICES 


FOR OPTIMAL NAVIGATION 

By John S. White 
Ames Research Center 


SUMMARY 


This paper presents the results of a study aimed at simplifying the 
equations used when applying the Kalman filter to space navigation. One fre- 
quently used method for generating the transition matrix is to integrate the 
perturbation equations of motion. An alternate approach considered here gen- 
erates approximate matrices directly from the gravitational potential. With 
this scheme^ the transition matrices can be generated much more quickly and 
with smaller computer storage requirements than before. Using these approxi- 
mate matrices for navigation is equally as good a technique as the original. 


INTRODUCTION 


In previous studies concerned with space navigation (refs. 1-6) the use of 
an optimal filter for processing the data has been considered. These studies 
have shown that such a filter permits accurate extraterrestrial navigation. 

Since the computations required are fairly complex^ it would be desirable to 
see what simplifications could be made without too much loss of accuracy. 

In this study we will consider the possibility of using simpler (and pos- 
sibly less accurate) methods for computing the transition matrix. This matrix 
is used in space navigation calculations to relate deviations of the trajectory 
at one time to those at another. Thus ^ 

x(t) = 9 (t;to)x(t^) 

where cp(tjtQ) is the transition matrix from tQ to t , and x(tQ) and x(t) 
represent the state at two different times. 

In the previous work the transition matrix has been obtained by integrat- 
ing the perturbation equations of motion (ref. 1^ appendixes B and E) . This 
approach utilizes a reference trajectory (obtained by integration of the true 
nonlinear equations of motion) to compute the coefficients of the perturbation 
equations and gives a transition matrix whose accuracy is limited only by the 
accuracy of integration. The computational procedure^ however^ is quite 
complex. 

Two possible means of simplifying this computation will be considered 
here. Firsts one can approximate the transition matrix by a power series in 
the two variables t and to^ and then use this series to compute the cp matrix; 
second^ one can derive a simplified expression for the matrix which is a 



fuaction of the partials of the gravitational field and the time increment- 
For both cases the necessary equations will be derived and results presented 
and compared with the results obtained by using the true matrices^ obtained by 
integrating the perturbation equation - 


MALYSIS 


In the previous studies (refs. 1-3) the transition matrix has been 
obtained by integration of the perturbation eq-oations of motion. The differ- 
ential equation governing the motions of the state variable x (composed of 
three position and three velocity components) is 

x(t) = F*(t)x(t) (l) 


The F* matrix can be partitioned into 3^3 submatrices ^ as 




F 


I 

0 


( 2 ) 


Where F is the gradient of the gravitational attraction. 


One of the properties of the transition matrix is that it satisfies the 
same differential equation as does x. Thus^ 

9(*^i^o) “ ^*("^ ) ) j T(*to^^o) ~ ^ ( 3 ) 


where double arguments are used to indicate the finishing and starting time. 
The solution of this expression can be expressed as a two-dimensional Taylor 
series in t^ and 6t^ where 5t = t - to with the origin for t and to 
assumed to be at the start of the trajectory. We then have for each element 
of the transition matrix 


CD (tjtp.) 



i=o j=o 


(^) 


The Cj_j are a set of appropriate coefficients that can be determined by 
a two-dimensional curve fitting process. That is^ one can evaluate *^o) 

along some specific reference trajectory by integrating the pert-urbation equa- 
tions and then evaluate a set of Cij that will cause the series to converge 
to the previously determined value of This process must be repeated for 

all elements of the matrix cp. 

This approach has two drawbacks. Firsts the shape of the cp surface is 
highly irregular as tp approaches the time of periapse. This is demonstrated 
in figures l(a) and (b)^ which are plots of two typical elements of cp. These 
are plotted for a constant time increment^ St^ and varying starting time^ to- 
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In contrasty plots of these elements versus St for fixed to (not sho*wn) are 
quite smooth. Thus^ for a good fit^ the series must contain high powers of 
to and therefore a large number of terms. Fitting the surface for only the 
midrange values of to one would require many fewer terms. Then ^ however, 
some other technique would be needed to obtain transition matrices at the end 
points of to* The other difficulty is that the are evaluated along a 

specific reference trajectory and would have to be re-evaluated for different 
trajectories . 


Development of F Series 


An alternate approach is to expand the solution of equation ( 3 ) in a one- 
dimensional Taylor series in St, giving 


<p(t;to) = cp(toJto) + q)(to;to)5t + cp(tojto) 


2 I 


(5) 


Equations for the derivatives of cp can be obtained appropriately differen- 
tiating equation ( 3 ) Q^ad evaluating the results at t = to • The derivatives 
of 9 are 

^o ) ^ ^ 

cp(to;to) = F* + F*F* ^ (6) 

= F* + 2F*F* + F*F* + 

etc . 

If the partitioned form of F* from equation (2) is substituted into 
equation (6) and these^ in turn, are substituted into equation (5)^ the follow- 
ing expression for cp results 


cpCtjto) = I + 


0 I 


F(to) 0 


St + 


F(to) 0 


j(to) F(t^^ 


2 I 


F(t ) 


F(t„) 


F(t^) + F(t_)^ 2F(t_) 


O'U 


5t^ 

31 


+ . 




^3 


(7) 
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•where 


cpi = I + F(to) 


5t^ 

2 ! 


+ 


j’(to) 



cp- = I 6 t + F(to) ^ + . . . 

^ 3 .' 

CP3 = F(to)6t + F(to) ^ + [F(to) + F(to)^] ^ 


cp^ = I + F(to) ^ + 2F(to) 


5t^ 

3 *’ 


where again 6t = t - t^. 


In this eq-uation^ F may he as complex or as simple as desired^ that is^ 

F may include only two-hody effects^ or may include additional terms due to 
ohlateness and other bodies of the solar system. Also it is apparent that one 
may include more or fewer terms in the series evaluation. For simplicity^ it 
would be desirable to use only the terms involving F^ dropping terms involving 
powers and derivatives of F. This gives an estimate of the transition matrix 
as 






F(t„)6t 


I + F(t^) 


21 


(8) 


Since F is a function of position in orbit,, this equation can be 
evaluated easily for any trajectory and any time along the trajectory. Thus,, 
no reference trajectory is required. Further^ it is relatively simple^ as 
discussed in the next section^ to determine how big 6t can be allowed to get. 


Development of ]Vlaximum Allowable Time Interval 

In using equation (8) there is some maximum value of bt,, beyond which 
the neglected terms in the series become important . This value of 6t can be 
determined by equating the first term dropped with the last term kept. Since 
the terms all involve matrices,, it seems reasonable to use the norm of the 
matrix in the equation. This leads to the following four equations. 
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l|F(to)|lSt - |l#(t^)|l ^ 

P(to)ll I? - 2ll^(to)|l ^ 

If these equations are solved for 6t^ one finds that the smallest 5t 
is obtained from the last equation^ and this value of 6t is used to determine 
a That is^ we wish to find a such that 

llF(to)ll = I llF(to)llSW (9) 

and then in evaluating equation (8)^ we will choose 6t < where k is 

less than unity- 

Let us define 


llFll = max 
1 ; J 


as being a reasonable norm to use in evaluating eq-uation (9)^ further,, let 

us restrict F to the two-body case- We then have (ref - 6^ eq. (6-51)) 


F = 




( 10 ) 


/ \ 

where (j. is the constant of gravitational attraction; r = (ri^ra^rsJ is the 
position vector of the vehicle with respect to the gravitational center; and 



In order to determine HfH^ we must choose that r which will maximize 
the element of F being considered. There are really only two types of terms 
to consider - diagonal and off diagonal. A typical diagonal term is 
fii = -(ii/r®)[l - (Sri^/r^)] which can he maximized if rx = r giving 
hlimaxl “ 2|x/r®. A typical off diagonal term is 

„ _ 3l^ rira 

^ r3 r^ 
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= (3/2) (n/rs) . Thus, 


which is a maximum when rx = giving |fiajj,axl 

choosing the larger of these tvo values ^ we have 


llFll = 2 

r3 


( 11 ) 


Having determined Hf 1| for use in equation ( 9 ) f 'we will now consider Hf| 1- 
Differentiating equation (lO) gives 


_ y 

3i^ 

r T 

rv 

^T 

-T- 

+ ^ ^ 1 


r 

r3 

rv 

rv 

rv 

^ yj 


( 12 ) 


where 


-T / \ dr 

V = (vx;.V2^V3) = 


dt 


and 


V = V 


The evaluation of 1 |fH is slightly more complicated since both r and v must 
be considered^ along with their scalar product. To simplify the evaluation of 
IIfII^ only values of f‘^v/rv = 0^ ±1 were considered^ and it was assumed that 
for other values of the product the maximum element would be nearly the same . 
Using the same general method as before^ we get 


l|F|l = 6 1 i 


( 13 ) 


Substituting equations ( 13 ) and (ll) into ( 9 ) and solving for &t 

max 

^max 2v 

Thus^ if 6t = r/2v^ the size of the terms neglected in equation (8) 
involving F) are comparable with those retained (involving F ) . It 
expected then that 6t will be chosen such that 6t < with k 

1.0 or less. 

A check was also made on the size of the second neglected term (not given 
in these equations)^ and it was found that for 6t = ^'tmax term at 

perigee is slightly less than the F term and becomes relatively smaller as 
the distance to the central body increases. Also_, since it involves an addi- 
tional power of 6t^ it decreases rapidly as 5t is reduced from 

In many cases ^ the desired &t (say from one observation to another or 
from the beginning to the first observation) is larger than thus 

two or more transition matrices must be combined to give the desired overall 
transition matrix. For example^ suppose we desire the matrix cp(ti5to) • 
However^ k^'tniax(‘^o) ^ ^"^a < “ '^o* readily generate where 

ta ^ to + from equation (8) . Then consider the time interval from tg^ to 

tx* Assuming that now k5tjj^2^(tg^) > tx - t^^ can obtain cp(*tajtg^)^ 
multiplying we get 


gives 

(1^) 

(those 

is 

being 
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9 (^ 1 ; to) = <p(tajta,)9(ta;to) (15) 

Of course^ this process imy be repeated as imny times as rec^uired to obtain the 
desired overall matrix. 


Effect of Two-Body Assumption 


The previous evaluation of F involved only two-body terms. One should 
also consider additional effects^ specifically^ earth oblateness and the per- 
turbing forces of the sun and moon^ and should determine the importance of such 
effects relative to the neglected terms in the power series expansion of F. 
Including these effects gives 


F = 


^^e __ 3rere^ 

re^ V re2 J 



3 


rs^ ) 


Ja- 


(I + 2M) { I 


5£e^ 
re^ J 


T - 
r- MTc 


I + 


T 

p Tprp 


2r-r ^ 


T 






where 


M = 


0 0 0 
0 0 0 
0 0 1 . 


(16) 


(17) 


It should be noted that this is the matrix form of the perturbation equations ^ 
which can be obtained by partial differentiation of the equations in appendix A 
of reference 1. We will consider the norms of each of the terms in (l6) 
individually and consider their relative effects . The first three terms have 
the same foim^ and the form of the norm has been already given in equation (ll). 
Thus^ 


Fo = 2 




"mil 


% 


IlF. 


= 2 ^ 
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The last term of eqviation (l6) reaches its maximum if r is all in the z 
direction vith the (3>3) term being the largest. This gives approximately 

llFjll = Ja2[3(-1>) - 5(1 + 2 - 7)J = 8 Ja^ 

These four components of ||fH are evaluated at four points along the earth-moon 
line to get representative values. The points used are at the earth's surface^ 
at 66^000 km from the earthy 66^000 km from the moon (the sphere of influence)^ 
and at the moon's surface. The distance from the sun was assumed constant at 
1 a.u.^ and the earth-moon distance was assumed to be 380^000 km. 

The results are tabulated in table I. As expected^ in the vicinity of the 
earth the Hf^H dominates^ while in the vicinity of the moon the ll^mll dominates . 
Thus^ if only a single two-body term is to be used for it must be switched 
between F^ and F^^ in the vicinity of the sphere of influence^ or alterna- 
tively^ both terms can be used at all times. The importance of these various 
terms will be discussed in the results section where ^ having determined a 
value for k^ we will be able to estimate the approximate size of the higher 
order terms of the series that have been rejected. 


DESCRIPTION OF SIMULATION 

In order to consider the effects of using approximate transition matrices^ 
a digital computer program was developed. This program is a modification of 
the one used by McLean (ref. 2). The principal modifications were made to the 
navigation system equations. These equations are discussed in the following 
subsection^ and are followed by a description of the trajectory used in the 
simulation. Other factors in the simulation are as used by McLean (ref. 2) and 
also by Smith (ref. 7)* 


Navigation System Equations 

When an approximate transition matrix is used in connection with the 
equations for the optimal filter^ the overall system becomes subopt imal. As a 
result^ the error covariance matrix obtained is only an approximation and is 
not a proper measure of the true estimation error. A true measure - that is^ 
a true covariance matrix of estimation error - requires additional equations. 
The optimal filter equations have been developed previously (refs. 1 and 2)^^ 
as well as the equations giving the statistics of subopt imal filters (ref. 7 ) • 
The necessary equations will be repeated here for convenience. It is assumed 
that observations are uncorrelated from one time to the next. 

^bviously^ the Kalman filter technique as applied to the navigation prob- 
lem is not optimal^ strictly speaking^ since the theory assumes a linear system 
and the navigation system equations are only approximately linear. However^ 
for the purposes of this report it is assumed that the linear approximation is 
exact and the covariance matrix of estimation error computed by means of the 
Kalman filter equation is^ therefore^ a true measure of system performance. 
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Consider first the true covariance mtrix of estimation errors. This 
matrix can he updated from the time of one observation to the time of the next 

where the - and 4 - superscripts refer to before and after an observation^ 
respectively; P'^(to) is given as an initial condition; and T(tj(;_j.2.^ t^;) is 
obtained by integrating the perturbation equations of motion. 

The estimated state x is also integrated from tj(; to tj^^-x* At the time 
of an observation the estimated trajectory is updated using an estimated 
weighting matrix as 


x"*” = X + Kgg-t(y - Hx ) (19) 

Here y represents the observation and H is the matrix of partial differen- 
tials relating deviations in the trajectory to deviations in the observation. 
Thus Hx~* represents the expected value of the observation based on previous 
data; is a weighting matrix which is expected to be an estimate of the 

optimum weighting matrix. 

As a result of this observation^ the true covariance matrix P is updated 
as 

P+ = (I - KestH)P~(l - KestH)^ + ( 20 ) 

This equation is valid for a nonoptimal system (see ref. 7 ^ (T^)) and gives 

the covariance matrix associated with the state estimation obtained in equa- 
tion (19) any K^st ’ 

Equations (18) and ( 20 ) could not be used in an actual system that used 
approximate transition matrices^ since the true transition matrices would not 
be available. For the purpose of this study^ however^ these equations were 
used to define the true state of affairs . 

The actual system will use the eq-uations for a true optimal system^ but 
will use an estimate of the transition matrix^ ‘^’est computed by means 
of the F-series (or any other desired method) . These equations will then pro- 
duce only an estimate of the covariance matrix and weighting matrix. 

The procedure starts by updating the estimated covariance matrix Pest 
by means of equation (18) ^ but using Test ^est-? where Pjst(-bo) = P+ho) . 
The estimated state x is updated to the appropriate time by integration and 
equation (I9) is used exactly as given to determine the new estimate. The 
weighting matrix required in equation (19) io update the estimated trajectory 
is given as 

Kest = PistH^(HPestH^ + Q)""” 
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is given as 


and Pgst 



^est ^est^est 


( 22 ) 


If the true covariance matrix P is used in equation (2l ) , the resulting 
Kest 'Will optimm. If this optimum Kgg^ is then substituted into equa- 
tion (20 ) , this equation will then simplify and have the same form as 
equation (22) . 

In order to get an indication of the nonoptimality caused by the use of 
the approximate transition matrix^ the true matrix^ is compared with the 
estimated matrix, ^est’ make this comparison, the quantities P^, and Py 
will be determined. These are defined as follows: The covariance matrix is 

partitioned into four 3x3 submatrices. The square root of the trace of the 
upper left submatrix is P^ and that of the lower right is Py. Thus, P^p 
represents the rms position deviation of the estimated trajectory with respect 
to the actual trajectory while Py represents the corresponding velocity 
deviations. Similar quantities from Pest-’ 'ttiat is, Prgst ^Vest-’ also 

determined . 


Trajectory 

In order to test the navigation accuracy using an approximate calculation 
of the transition matrix, a sample circumlmar mission was chosen for a refer- 
ence. The particular reference trajectory chosen has a 7O.68 hour flight time 
to the moon with a perilune altitude of I87 km. It has a free return to the 
earth arriving at perigee at l44.87 hours after launch and entering the earth’s 
atmosphere near the center of the entry corridor. The navigational observa- 
tions used were theodolite measurements, which determine the right ascension 
and declination of the earth or moon as seen from the vehicle. Appropriately 
chosen sextant data would give essentially the same results. There were 
45 observations made on the outbound leg of the mission and 39 "the return 
leg. The observations were grouped in a manner very similar to that shown in 
figinre 2 of reference 7; that is, observations of the earth were grouped 
together, as were those of the moon, and some of both types were made shortly 
before each of the three expected velocity corrections on each leg. 


RESULTS AM) DISCUSSION 


The purpose of this section is to discuss the operation of the navigation 
system using approximate transition matrices and to consider the effects of 
the various simplifications. First of all, we will discuss the accuracy of the 
F-series transition matrices. Then we will consider the step size required for 
use with the F-series expansions, and will compare the resultant covariance 
matrices with those obtained using the exact transition matrix in a truly 
optimimi system. The effect of gravitational perturbations is then considered 
followed by a brief discussion of computer storage and time requirements. 
Finally, the effect on guidance will be considered. 
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Transition Matrices 

In the computer program the true transition matrix^ calculated hy 
integrating the perturbation equations^ -was available^ as well as the F-series 
approximation. These two matrices were compared (using the norms of the sub- 
matrices)^ and the percentage difference between the estimate and the true 
transition matrix varied from 1 to 0.1 percent for long inteiTvals where many 
transition matrices must be multiplied together and from 0.1 to 0.001 percent 
for short intervals where a single matrix can be used. McGee (ref. 8) states 
that adequate navigation can be performed if the transition matrices for short 
intervals are calculated using only l6 bits or O.OO3 percent error- The 
F-series transition matrix has about this error for the short intervals and 
thus one would expect that it would provide adequate navigation^ as indeed it 
does . 


Effect of Step Size 

In order to determine the effect of using various time intervals in calcu- 
lating the transition matrix from the F-series^ a comparison was made between 
the approximate covariance matrix of errors in the estimated trajectory^ 
and the true covariance matrix. The F-series itself was calculated omit- 
ting the sun and oblateness terms^ but including at all times both the earth and 
moon terms to provide a smooth transition from earth-centered to moon-centered 
conics . 

A study of Pp and P^ from both the true and estimated covariance matri- 
ces will then show the adeqmcy of the transition matrix approximation. A 
plot of ^rest/^^ both outbound and return legs of the cir- 

cumlunar trajectory is shown in figure 2 for four different values of k. 

Values of P and Pest available only after an observation. The gaps in 

the curves represent periods of no observations. Ideally^ of course^ these 
ratios shoTiId be constant at unity. 

It can be seen in figures 2 (a) and (b) that, at least until perilime, a 
value of k = 0.1 gives quite good results and that the accuracy falls off as 
k increases, so that at k = 1.0 the results are poor. After perilune, how- 
ever, all of these curves deviate from unity, and data not shown indicate that 
this trend continues so that at perigee the ratio is quite far from unity. 
However, if the run is restarted at perilune with P = Pest nominal values, 
then as shown in figures 2(c) and (d) , the results are very similar to those 
of the outbomd leg, with k = 0.1 quite good and k = 1.0 fairly poor. 

It is apparent, however, that smaller values of k are required for 
adequate results after passing perilune. If the F-series is calculated at 
every integration step, which implies that k is small, the ratio of Pest/^ 
remains close to 1.0 for the complete round trip. Thus, these smaller values 
of k in the vicinity of the moon allow a successful perilune passage. How- 
ever, this gives a considerable increase in the amount of computation. If the 
vehicle is returning from a lunar orbit, it would be natural to restart the 
Pest matrix at its nominal value so that extremely small values of k would 
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not be required. In a fly-by mode^ however^ this would not be natural. 
Nevertheless^ at perilme a nominal P matrix could be inserted in the com- 
puter^ which would allow successfilL computations during the retiorn flight. 

These results indicate that the assumptions made in determining dt^^^y 
(eq. (l^)) are reasonable^ and that a value of k = 0.25 gives reasonably 
accurate results. This value is used as a standard in the remainder of the 
discussion. 

Using this value of k^ one must multiply several short transition matri- 
ces together to cover some of the larger required intervals. For instance^ 

30 short term matrices were required to generate the first matrix from t = 0 
to 2.5 hours ^ and 15 were required during the last 48 hours to perilune . Most 
of the time a single matrix was adequate to cover the time interval between 
observations . 


Comparison With Optimum System 

It is of interest to compare the results of using the approximate transi- 
tion matrices with the results of using the true matrices in a trixly optimum 
system. This comparison is made in figures 3 and 4. In figure 3 ^ the true 
transition matrices and the optimal K matrices were used throughout. Thus 
Pegt ^ equal and optimum. The curves represent P^^ and and the 

points represent the navigation error of the estimated trajectory for a single 
typical flight. The navigation error is supposed to be a typical member of 
the ensemble for which P is the rms value . It can be seen that the points 
could easily represent a typical member. 

In figure 4 the estimated transition matrix was used ^ and Prest ^"^est 
are plotted along with the associated navigation error. The curves in fig- 
ures 3 and 4 are similar with only small differences^ indicating that Pest 
is a good measure of the optimum P. The navigation error points are also 
similar. This is to be expected since identical sequences of simulated obser- 
vation errors were used in each case^ and the covariance matrices are also 
similar . 


Effect of Two-Body Calculations 

In the previous results^ both the moon and earth gravitation terms were 
used at all times to compute the F-series. It might be simpler to use only 
one of these terms and to switch terms at the sphere of influence . These two 
techniques are compared for k = 0.25 in. figure 5* The plots of Pest/^ 
the two cases are essentially identical from launch to the sphere of influence 
of the moon. Following this,, there is some difference^ especially in fig- 
ure 5(b),, which tends to decrease with time. This same effect is also notice- 
able during the return trajectory. In general,, either technique will work 
adequately,, but the three-body technique gives slightly better results. 
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Effect of Otlateness and Solar Gravitation 

It is also of interest to ascertain the solar and ohlateness effects. In 
table I the size of these terms is given along with those of the earth and 
moon. The sim's effect is always less than that of either the earth or the 
moon^ and so will have considerably less influence than that shown in figure 5 
for the difference between two- and three-body calculations. The oblateness 
term is largest at perigee, at this point, for k = 0 . 25 , the second term in the 
series, which is being neglected, is about 0.6(l0)^, or about 4o times larger 
than the oblateness term. Since the relative size of the oblateness term 
decreases rapidly on leaving perigee, it is apparent that it will have less 
effect than truncating the series. In conclusion, it appears that the solar 
gravitation and earth's oblateness will have very little effect on the overall 
results . 


Computer Storage and Time Comparison 

Since the F- series calculation of the transition matrix appears to give 
reasonable results , it is of interest to determine the computer storage and 
computational time requirements for comparison with the corresponding require- 
ments of the technique where the pertiorbation equations are directly integrated. 
Computer programs are available for both techniques, and a limited amount of 
data is available concerning storage and timing. However, both programs are 
written in Fortran IV, and no optimization of either has been attempted. Thus, 
the storage and time requirements given here must be considered approximate . 

For each of these two programs the storage specifically required for the cal- 
culation of the transition matrix and the time required for performing these 
calculations during the outbound leg of the mission were determined and are 
presented in the following table: 

Storage , words Time , min 

F-series 600 (k = O.l) l/k 

(k = 0.25) 1/5 

Pertiirbation integration l400 3/4 

It can be seen that the F-series takes somewhat less than half the 
storage required for the perturbation integration technique. The time required 
to integrate the perturbation equation is constant at about 3 /'^ of ^ minute; 
whereas the time required for the F-series calculations is a function of k. 

In the extreme where the F-series is evaluated at every integration step, the 
time required is about 3 /^ minute. Thus, the F-series approach results in 
savings in both storage and time. 


Effect on Vehicle Guidance 

It has been shown that an approximate transition matrix can be used 
satisfactorily for navigation in cislunar space, that is, for determining the 
vehicle's present position and velocity. Once the navigation is completed, the 
results are used in one of several available guidance procedures to determine 
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wiiat, if any^ corrective maneuver is required so that the vehicle will 
satisfactorily complete the mission. Using the transition matrix for guidance 
as well as navigation has heen st-udied in previous research (refs. 2^ 6). In 
those studies it was ass-umed that the vehicle woilLd be required to arrive at a 
desired specified point at a particular time. The transition matrix^ from the 
present to the end pointy was obtained (ref. 2) by miiltiplying many short 
individml matrices together. This guidance scheme was used in the present 
study in two ways: Firsts the approximate transition matrices were used for 

navigation only^ and the true transition matrices were used in the guidance 
scheme. In this case^ the standard deviation of the miss at arrival was about 
the same as previously reported^ namely^ about 2 km in altitude and 20 km in 
range. G?his indicates that the approximate transition matrix provides adequate 
navigation. Second^ it was desired to see the effect of using the approximate 
transition matrices in the guidance section also. In a study of word length 
requirements (ref. 8)^ it was concluded that this guidance scheme would be 
usable with only l6~bit accuracy in the computation of the transition matrices. 
The approximate matrices used here have about l6~bit accuracy for the majority 
of the time^ but at the beginning and at the end^ the accuracy decreases to 
about 6 bits. When these approximate transition matrices were used in the 
guidance scheme^ the standard deviation of the altitude miss was about 1000 km^ 
indicating that these transition matrices are unsatisfactory for this purpose. 

Thus^ one can conclude that these approximate transition matrices are 
quite satisfactory for use in navigation^ but that a guidance scheme should be 
selected which does not use them. 


CONCLUSIONS 


This report shows that a simplified technique for computing the transition 
matrices is available which provides adequate navigation in cislmar space and 
which uses less computer storage and less computation time. These transition 
matrices are computed as a function of the gravitational attraction along the 
estimated trajectory and so do not require a reference trajectory. 

The gravitational attraction can be computed using the standard two-body 
equations. The effect of the earth’s oblateness and the sun can be neglected. 
In going to the moon^ the moon’s attraction can either be included all the 
time^ or the gravitational center can be switched at the lunar sphere of 
influence . 

The maximum time over which a single calculation of the gravitational 
attraction can be used to generate the transition matrices is easily calculated. 
Where a longer time interval is required^ several individual transition 
matrices can be computed and then multiplied together to obtain the overall 
desired matrix. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Fields Calif. ^ Feb. 28^ I 966 
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APPENDIX 


DEFINITION OF SYMBOLS 


F 

F* 

llFll 

f 


IJ 


H 

I 

J 

K 

k 

0 

P 


- V 


Q 

r 

r 

t 

V 


equatorial radius of tlie earth. 

coefficient of the ijth term of a two-dimensional power series 
gradient of the gravitational attraction vector 

, the matrix of coefficients in the differential equation of 


0 I 
F 0_ 

" motion of the state vector 


norm of the F matrix 

the ij element of the F matrix 

matrix of partial derivatives of the observed quantity with respect to 
the state variables 

identity matrix of suitable dimension 

coefficient of the second harmonic of the earth potential 
optimal weighting matrix 

constant used in determining allowable step sizes 
null matrix of suitable dimension 
covariance matrix of estimation errors 

rms position deviation of the estimated trajectory with respect to the 
actual trajectory 

rms velocity deviation of the estimated trajectory with respect to the 
actual trajectory 

covariance matrix of observation errors 

Irl 

T 

(ri^r^^rs) ^ position vector of the vehicle 
time 

T 

(vi^va^vs) ^ velocity vector of the vehicle 
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X 


6t 

U 

9 


9i^92^ 

93.94 


O 
( )" 

( )-^ 


e 

est 

i.j A 

J 

m 

o 

s 


+ 


state vector representing position and velocity deviations from the 
reference 

state vector of estimated deviation from the reference 
time difference over which the transition matrix is used 
maximum allowable value for St 
constant of gravitational attraction 

transition matrix relating deviations at one time with those at 
another 

submatrices of cp 


Rotational Conventions 
time derivative of ( ) 
transpose of matrix ( ) 
inverse of matrix ( ) 

Subscripts 

earth 

estimated 

indices 

earth oblateness 
moon 

initial condition 
sun 

Superscripts 

value before an observation 
value after an observation 
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TABLE I.- TEE NORM OF THE COMPONENTS OF F VERSUS LOCATION 

IN CISLUNAR SPACE 


Vehicle location 

llFell 

IlFmll 

IlFsll 

iifjii 

Perigee 

66^000 3 an from earth 
Lunar sphere of influence 
Perilune 

2.44(10)^ 

2.20 

•193 

.115 

l.4i(io)“® 
2. 36 (10) -3 
.270 

1.49(10)^ 

6.21(10)“^ 
6. 21(10) 

6. 21(10) 
6.21(10)-^ 

1.57(10)2 
1.52(10)“3 
4. 65 (10) -7 
1. 98(10) 
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(a) Typical diagonal element. 

Figure 1.- Typical elements of the transition matrix versus initial time with 

final time 1*2 hoixrs later. 





(b) Typical off diagonal element. 


Figure 1.- Concluded. 
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(a) Position^ outbound, leg. 
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Figure 2.- Comparison of estimated, and. actual covariance matrices for various values of k. 
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(Td) Velocity,, outlound leg. 



Figure 2.- Continued. 
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(c) Position^ return leg. 
Figure 2.- Continued. 
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(d) Velocity^ return leg. 


Figure 2.- Concluded. 
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Figure 3-“ Navigation error and its standard deviation for the optimum, case 
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(a) Position. 

Figinre 4.- Navigation error and its standard deviation using approximate 

transition matrices^ k = 0 . 25 • 


27 



standard deviation Pw km/sec 


2.0X10""^ 



0 10 20 30 40 50 

t, hr 


(b) Velocity. 
Figure 4.- Concluded. 
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conducted so as to contribute ... to the expat2sion of human knowl- 
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shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof.” 
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